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We study the effects of noise on a recently discovered form of intermittency, referred to as in- 
out intermittency. This type of intermittency, which reduces to on-off in systems with a skew 
product structure, has been found in the dynamics of maps, ODE and PDE simulations that have 
symmetries. It shows itself in the form of trajectories that spend a long time near a symmetric 
state interspersed with short bursts away from symmetry. In contrast to on-off intermittency, there 
are clearly distinct mechanisms of approach towards and away from the symmetric state, and this 
needs to be taken into account in order to properly model the long time statistics. We do this by 
using a diffusion-type equation with delay integral boundary condition. This model is validated by 
considering the statistics of a two-dimensional map with and without the addition of noise. 
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I. INTRODUCTION 

Many dynamical systems of interest possess symme- 
tries which force the invariance of certain subspaces. A 
great deal of effort has recently gone into the study of 
such systems, in particular studying the behaviour of the 
attractors near their invariant subspaces on varying a pa- 
rameter §, §, §, |, §, §, @, || . This has included the study 
of systems with both normal and non-normal parameters 

Such systems show a variety of novel phenomena in 
their dynamics. In particular, systems with normal pa- 
rameters that are of skew product type (namely those 
where the transverse dynamics does not affect the dy- 
namics tangential to the subspace) can show on-off in- 
termittency, which occurs as the result of the transversal 
instability of an attractor, usually chaotic, in the invari- 
ant subspace whose trajectories get arbitrarily close to 
the invariant subspace, while making occasional large de- 
viations away from it || 0, S. On-off intermittency can 
be modelled by a biased random walk of the logarithmic 
distance from the invariant subspace || [|, || . 

On the other hand, systems with non-normal parame- 
ters which do not have skew product structure can show 
other dynamical phenomena in addition to those present 
in skew product systems. These include a type of inter- 
mittency referred to as in-out intermittency p0|| ; sim- 
ilar effects were noticed independently in a number of 
models jll], Q. Examples have been recently found in 
PDE models of surface waves |13 and in a problem of 



chaos control in the confinement of magnetic field lines 
in toroidal fusion chambers [fl4| . In the original formula- 
tion of in-out intermittency dropping the condition that 
a chaotic attractor is necessary in the invariant subspace 
turned out to be an important ingredient [ fTo] , and this 
has subsequently been shown to lead to further novel 
phenomena JlEfl . 

This type of intermittency is best characterized by con- 
trasting it with on-off intermittency. Briefly, let Mj be 
the invariant subspace and A the attractor which exhibits 



either on-off or in-out intermittency. If the intersection 
Ao = A n M] is a minimal attractor then we have on- 
off intermittency, whereas (in the more general case) if 
Aq is not necessarily a minimal attractor, then we have 
in-out intermittency. In the latter case there can be dif- 
ferent isolated invariant sets in Aq associated with at- 
traction and repulsion transverse to Aq, hence the name 
'in-out'. Another difference is that, as opposed to on- 
off intermittency, in the case of in-out intermittency the 
minimal attractors in the invariant subspaces do not nec- 
essarily need to be chaotic and hence the trajectories can 
instead shadow a periodic orbit in their 'out' phases ]To|| . 
A schematic representation of this scenario is depicted in 
Figure [l[ 

There is now a good understanding of the statistical 
properties of on-off intermittency [l7| and some 

properties of in-out Jiof] intermittency where they differ. 
Numerical support has also been obtained for both on- 
off |, H § ( for experimental evidence see [^7], |l) ) and 



in-out ||24], B5J . 

Both of these types of intermittency rely on the pres- 
ence of invariant subspaces. In real systems, however, 
invariant subspaces are only expected to occur approxi- 
mately; either as a result of the lack of precise symmetry 
or due to the presence of noise. This has motivated a 
number of studies of the effects of noise on the statistics 
of on-off intermittency g|, H 13= H EH !§■ 0ur aim 
here is to make an analogous study of the effects of noise 
in the case of in-out intermittency, by making a contin- 
uum version of the Markov model considered in fltj|| in 
order to highlight the similarities and differences. We 
do this by considering an analogue of the drift-diffusion 
model employed by ^, [l?], [l^, ^| for on-off intermit- 
tency. 

The structure of the paper is as follows. In Section || 
we derive and analyze a model of in-out intermittency 
that consists of a drift-diffusion equation with delay in- 
tegral boundary conditions, based on extracting the im- 
portant information from a dynamical model. We also 
discuss how to estimate parameters in the drift-diffusion 



2 




if it has a basin with positive Lebesgue measure, such 
that any smaller invariant set has a basin with smaller 
measure. An attractor is minimal if it contains no proper 
subsets that are attractors). 

Suppose now that j(t) is a typical trajectory in the 
basin of A, such that the w-limit set of j(t) is the attrac- 
tor A. We assume that A° C Aq is a Milnor attractor 
contained within A for /|mj- We assume also that the 
only transversely stable set in Aq is some A 1 C Aq which 
is a repeller for /|mj- Each of the invariant sets A l '° is 
assumed to support invariant measures \x l and \x° that 
govern the behaviour of typical trajectories in A on the 
approach to A l '°. 

We assume that A 1 is transversely attracting on aver- 
age (i.e. its largest transverse Lyapunov exponent (L.E.) 
with respect to fi l is X l T < 0) and A° is repelling on aver- 
age (i.e. its largest transverse L.E. with respect to n° is 
positive). We refer to A 1 as the 'in' dynamics and A° as 
the 'out' dynamics for obvious reasons; see Figure [l] for 
a schematic representation. 



FIG. 1: Schematic diagram showing a typical trajectory that 
is in-out intermittent to an invariant submanifold Mr, show- 
ing the 'in' and "out" phases. The invariant submanifold Mr 
contains an invariant set that decomposes into a transversely 
attracting chaotic saddle A 1 and a transversely unstable peri- 
odic orbit A" that is an attractor within Mi. The 'out' phase 
is defined by the trajectory being within the isolating neigh- 
bourhood U° of A° . In this case, the 'in' phase is modelled by 
a random walk in the logarithmic distance from Mi whereas 
the "out" phase shows uniform exponential growth away from 
Mi, shadowing the unstable manifold of A". 



model. Section III adapts these to include additive noise 
in the transverse variable as well as in the tangential 
variable. We predict new transitions in the dynamics on 
adding noise to the tangential variables. Section IV dis- 



cusses the estimation of the parameters in the model and 
obtains scalings and transitions on changing the noise 
amplitude. These predictions are tested on a planar map- 
ping given in (To) . Finally Section ^ gives a discussion 
and interpretation of the results. 



1. Identifying the 'in' and 'out' phases 

Suppose now that we have a projection II : R" — > Mi 
and a neighbourhood U° C Mj containing A°, such that 
f (U°) D U° (i.e. it is absorbing for f\ Ml ) ||. We iden- 
tify a point x in the phase space as being on 'out' phase 
if II(x) <E U° and as being on the 'in' phase otherwise. 
Note that there is an arbitrary choice of neighbourhood 
U° and projection II; however we will be interested in 
statistical properties of the 'out' phases that are inde- 
pendent of these. 

For concreteness (and to correspond with examples 
studied later) we assume that the dynamics on A° is pe- 
riodic and the dynamics on A 1 is chaotic (with many 
ergodic measures supported on A 1 ). However, in prin- 
ciple the same type of model applies as long as at least 
one of A 1 or A° is chaotic. We are interested in modelling 
the asymptotic fluctuation of the distance of some typical 
trajectory j(t) from Mj. Suppose we have a function 



II. MODELLING IN-OUT INTERMITTENCY 

A. The dynamical model of in— out intermittency 

Suppose that we have a dynamical system that evolves 
on R", such that some subspace Mj (dim(M/) < n) of 
R n is dynamically invariant. For definiteness we consider 
a dynamical system generated by iterating some smooth 
map / : R" -> R", in which case f(Mj) CM/. If there 
is a minimal Milnor attractor A for this system such that 
A n Mi = Aq is not a minimal Milnor attractor for the 
system restricted to Mj, then we say the attractor is in- 
out intermittent jl0| . (Recall that A is a Milnor attractor 



y(t) = nj(t) (1) 

where it : R" — > R is a smooth function such that 
7r _1 (0) = Mi. Then we say ir projects the phase 
space onto the transverse variable y. Clearly we have 
liminf^oo \y(t)\ = but limsup t \y(t)\ > 0. Moreover, 
such a transverse variable will, because of invariance of 
Mi, spend arbitrarily long times near y — 0; the so- 
called laminar phases. The object of this paper is to give 
a statistical description of the behaviour of a generic y(t) 
measuring the distance from M/ for in-out intermittent 
dynamics. 
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B. A Fokker-Planck model for in— out intermittency 

1. In terms of a logarithmic transverse variable, z. 

We start with a transverse variable y = tt(j) and set 
z = — In we model the behaviour of z as follows. Dur- 
ing the 'in' phase, we model the behaviour as though it is 
a linear skew product forced by the chaotic 'in' dynamics 
and we assume, by an appropriate scaling, that \y\ < 1 
for all time. We model the behaviour as a drift-diffusion 
process in z > with drift — > per unit time and 
diffusion f3 2 per unit time subject to reflection boundary 
conditions at z = 0. We assume that the trajectory leaks 
onto the 'out' phase at a rate e > per unit time (this is 
given by the most positive tangential L.E. of /x.;). 

On the 'out' phases we assume that there is a fixed 
linear expansion forced by the periodic 'out' dynamics. 
This translates to a deterministic growth in the z variable 
at a rate — < 0. Once z reaches we assume that the 
trajectory is forced to re-inject to the 'in' dynamics. For 
convenience, from here on we define 



c 
A 



A 
A' 



at z = 0, which corresponds to reinjecting trajectories 
reaching z = of the 'out' chain back into the 'in' chain. 
If we define the total overall probability of being in the 
'in' and 'out' chains by 

= &(t) + $°(t) (8) 

then we have 

d<S> , , , , 

(9) 



dt 



00 'dP dR\ 

-z- dz 



dt dt 



■ 'dT OR n . , 
eP + c— + eP\ dz 



io \dz dz 
= T(0,t) + cR(0,t) 
= 0, 

implying that $(i) is a constant. We therefore stipulate 
that by normalization 

*(*) = i (io) 

for all t. Thus the Fokker-Planck model of in-out in- 
termittency (in the absence of noise) is the closed sys- 
tem consisting of the linear equation (|^) for P(z, t) on 
z € [0, oo) subject to the delay integral boundary condi- 
tion (|7T) and normalization condition (fol). 



and note that c > 0, whereas A can be positive or negative 
in what follows. 

Let the probability density at time t of the distribu- 
tion of z values in (z, z + dz) on the 'in' phase be given 
by P(z,t)dz and those on the 'out' phase be given by 
R{z, t)dz. Our model translates to a forward Kolmogorov 
equation for P of the form 

dP 
~di~ 



dr 

dz 



(2) 



where the second term on the right hand side represents 
the leakage into the 'out' phase and 



dP 

T(z,t) = f3 2 /2?- + \P 

dz 



(3) 



represents the flux of trajectories at (z, t). The dynamics 
for R on the 'out' phases is simply given by the hyperbolic 
equation 



— - c— r 

dt dz , 



eP. 



This equation can be solved exactly to give 

R( z , t) = - f P[X,t- — - ) dr. 
c.L \ c 



(4) 



(5) 



which is unique up to the addition of an arbitrary func- 
tion £(x + ct). The total probability of being in the 'in' 
or 'out' chain is then given by 

/•oo />oo 

® l (t)= P(z,t)dz, $°(<)=/ R{z,t)dz(6) 
Jo Jo 

respectively. We assume also that 
r(o, t) + cR(Q, t) = 



(7) 



2. In terms of the transverse variable, y. 

The drift-diffusion model in the variable z can be trans- 
lated into one for the original transverse variable y = er z 
as follows: Let the probability density at time t of the 
distribution of y values in (y, y + dy) on the 'in' phase be 
given by Q(y,t)dy and those on the 'out' phase be given 
by S(y, t)dy. Note that (assuming y > 0) 



P = 



dz 



Q = yQ, R = 



djl 

dz 



S = yS 



(11) 



and so the system governing Q(y,t) and S(y,t) is given 
by 



s = — 

cy 



w,t + 



In y — In w 



dw, 



with the boundary conditions given by 



1 



[Q(y,t) + S(y,t)] dy 



io 

= ^-(y-i-iyQ)) ~ x vQ + s at v 

2 dy 







Observe that in terms of these variables we have 
$'= / Q(y,t)dy, / S(y,t)dy. 



Moreover we find 



(12) 
(13) 

(14) 
(15) 

(16) 

(17) 
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C. Stationary distributions; noise free 



Steady solutions P(z, t) — p(z) of (0) will satisfy 



ep : 



(18) 



with boundary conditions given by (§) and (JToj) . This 
can easily be solved to give a solution 



p(z) = Ae^- Z + A + e^+ 2 



where 



-X± VA 2 + 2e/3 2 
J 2 



(19) 



(20) 



Note that if e > 0, then a/A 2 + 2e/3 2 > |A| and so /i_ < 
and n+ > always, as long as e > 0. A similar result 
was found for the Markov model of in-out intermittency 
discussed in [fbl. We therefore write 



M = M- 



(21) 



and note that the only solutions of (|18|) with finite mass 
are such that A + = 0. Calculating the stationary mass 
in the 'out' chain we have 



r(z) = Be» z 
where —c\iB = eA and so B 

A* ' A* 

Normalizing so that the total mass is unity we have 



(22) 

A. This means that 
(23) 



1 = ^ + 5 
A* A* 



which gives 



C[l 



B = 



Cfl 



CfJ, 



C[l 



(24) 



(25) 



thus ensuring that the boundary condition is also satis- 
fied 



r(o, t) + cR(o, t) 



Ae^ = 0. 



(26) 



These steady exponential distributions correspond to al- 
gebraic distributions for the Q(y) and S(y). In Section 
1VA we discuss how the free parameters in this model 



can be estimated from the dynamical data. 



D. Contrasts with on— off intermittency 

Note that one could take a simpler dynamical model in 
the form of a drift-diffusion equation but with no differ- 
entiation into 'in' and 'out' phases. This is equivalent to 
assuming that e — in the model (i.e. reducing it to an 



on-off process) and leads to an exponential probability 
distribution of the form 



p(z) =4e ? 



(27) 



However, estimation of the constants A, (3 presents a 



problem as we will discuss in Section IV B| . 

For in-out intermittency, the 'out' phase is distinct 
from other invariant sets that we may choose within the 
invariant subspace in the following sense. There are con- 
stants E,F > such that if any trajectory enters the 
'out' phase at a distance d from the invariant subspace 
then there is a minimum residence time in the 'out' phase 
given by 



.Find. 



(28) 



In particular the minimum 'out' phase residence time 
goes to infinity as d — > 0. 



III. A MODEL FOR IN-OUT INTERMITTENCY 
WITH ADDITIVE NOISE 

The model of Section |l| can now be easily generalized 
to model in-out intermittency in the presence of unbi- 
ased additive noise. At first we will investigate the case 
where the noise is added only to the transverse variables 
and later examine the case of noise added to tangential 
variables. It appears that noise in the transverse vari- 
ables affects scalings in a regular manner; whereas noise 
in the tangential variables can lead to transitions as 'in' 
and 'out' phases merge. 



A. Additive noise in the transverse variables 

In presence of additive noise in the transverse variables, 
we use a similar approach to that of Ashwin & Stone |lj] 
and Venkataramani et al. jl7j to obtain a Fokker-Planck 
model of the form 



dQ 
dt 



-A 



d(Qy) 
dy 
a 2 d 2 Q 
2 dy 2 



0^d_ 
2 dy 



d(Qy) 
dy 



(29) 



which is similar to (|12|) apart from the diffusion term at 
a rate tr 2 /2 corresponding to the additive noise. 



1. Steady state with additive noise 

We obtain a steady state probability density distribu- 
tion in the in-out case, by calculating the contributions 
from both 'in' and 'out' chains separately. 

For the contribution from the 'in' chain we proceed by 
looking at the steady state counterpart of the equation 
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( p9| ) obtained by demanding Q(y,t) = q(y), 

d(yq) 8 2 d ( d{yq)\ 
- X ^y- + Ydy-{ y ^y-) + 



(30) 



a 2 d 2 q 



eq = 



2 dy 2 

which can be written in the form 

^ i . > n . o i — % A- ( — B 2 — \ \ r 

To solve this equation we recall that the case with e = 0, 
corresponding to on-off intermittency, is solvable explic- 
itly (see [[I?]) with the solution 



q(y) = A(^y 2 +a 2 )^ 



(32) 



where £ = In the case of in-out intermittency with 
e / 0, we proceed by employing the (singular at a = 0) 
change of variable 



to rewrite the equation (^l]) in the form 



(33) 



/3 2 - 2A - 2e 
I/3 2 



g = 0. 



This equation can be solved in terms of hypergeometric 
functions 



q — A{(3,a,\, e)F( a+j0l _ tl / 2}X ) 



where 



a± 



1 - A ± VA 2 + 2B 2 e _ 1 + p± 



2/3 2 

with n± as before, and 



(35) 



(36) 



Ft 



(a,b,c,z) 



r(c) r^(i-t) 



c-b-l 



T(b)T(c-b)J {l-tz) a 
o 



-di. (37) 



If |i| > 1, x < and i?ea;(a - b) = 2^A 2 + 2^ 2 e > 0, 
this solution can be approximated by using 



F 



(a + ,a-,l/2,x) 



U-^J x r(a+)r(i/2-Q_) K ' 

K(l-x)- a -, 



for some constant K. This then gives 
g(y)«A( CT 2 +/3V)" 



(39) 



where A = A(B, a, A, e) is the normalisation constant, 
and n — /i_ as before, is valid for small \y\. 



For the 'out' chain, the steady state probability density 
distribution can be obtained from (|l^) by employing the 
steady state form of Q in the integral in this expression 
to compute s[y) 

(40) 



s(y) = — / q(w)dw. 

cy Jo 



The overall steady state probability density distribution 
is then given by the sum of these two contributions q(y) + 

s (y)- 



2. Scaling of the mean first crossing time 

To begin with, we recall that for unbiased noise on the 
variable y, the mean is clearly 



E(y) = 0. 



(41) 



To determine the mean crossing time through y = 0, 
we require A which can be computed assuming there are 
no- flux boundary conditions at y = ±1 (with linear be- 
haviour up to this point), to give 



P(y)dy=l. 



(42) 



\v\<i 



Away from a blowout bifurcation point, it is not so 
easy to obtain an explicit expression for A and thus for 
the variance var(y). However, in the low noise limit, 
a <C 8, we can approximate the stationary density (|39| ) 
by the continuous function 

J Act- 1 -^ for \y\ < a/(3, 



A{py)- 1 -" for */0<\y\<l, 
and s(y) is calculated from ( [40| ) as 



s(y) S3 



-i-n 



iPv) 



/3y 



for 1 2/| < er//3, (- 44 ) 
for \y\ > a//3. 



In order to compare these results with those in the case of 
on-off, we proceed by computing the value of the normal- 
isation constant A. This is given by computing ( [42] ) with 
P{y) — q{y) + s(y) and gives, for this approximation, 



A = c^ 2 I3 [(e - c/z)(/T" - a~") 



(45) 



In the limit e 

case 

A = 



-<7-"/i(elng(l+/*) + (c+e)/i)) . 

this expression reduces to the on-off 



ftp 



o-"(l + //) - 



(46) 



equivalent to that found in Ea |. Using the approxima- 
tion that the stationary distributions g(y) and s(y) are 
approximately constant for |t/| < er, the instantaneous 
flux from y > to y < can then be estimated as 



(47) 
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where, as in Jjj|, we have assumed that with unbiased 
additive noise approximately half of all initial points in 
[0, a] will cross over within the next time unit. We can 
compute J 7 , by using the above solutions (|||) and (|44| ) 
together with (f45[), to obtain 



T = a-^^ 2 (3(c + e) [2((e - c/x)(/?-" - o-") 
-<r-MeInf(l + /i) + (c + eW) 



(48) 



which in the limit e — > reduces to the on-off formula 



(49) 



This allows the calculation, in the in-out case, of the 
mean first crossing time given by M = 1 j T analogous to 
M- i.e. 



M = 



(( e -c M )((|)-P-l) 



- eAi ln|(l + /i) + (c + e) M 2 )) 



(50) 



Considering the case where e is asymptotically small (i.e. 
small leakage of the 'in' dynamics) we have 



(3 




0(e,p) 



(51) 



Furthermore, if we are close to marginal stability on the 
'in' chain, A = 0, we have p, = —^if- and so in the limit 
£->0we recover the expression in |18| , namely 



M 



(3 (ip 




- 1 



(52) 



with an order y/e correction. Similarly, expressions can 
be obtained for other limiting cases; we give one such 
scaling with the numerical results in a later section. 

We have attempted to find the scaling of the mean 
laminar length with noise intensity, analogous to pc| for 
on-off intermittency. However, the need to distinguish 
between the dynamics of the 'in' and and 'out' phases 
means that we cannot easily reduce the problem to a sin- 
gle ordinary differential equation with the consequence 
that we have so far not been able to obtain an expres- 
sion as compact as that for on-off. However, in principle 
the Fokker-Planck model (|2^) contains all the necessary 
information to compute this. 



B. Added noise in tangential variables 

It has been noted that the addition of noise to 
tangential variables in the case of on-off intermittency 
has only a minor effect on the dynamics. This can be 
understood if the attractor within Mi is stochastically 
stable, i.e. if the probability density with noise limits to 
the probability density of the natural measure in the case 
of no noise. In the case of in-out intermittency, on the 



other hand, there will be a threshold of noise amplitude 
beyond which the fine structure in the invariant subspace 
is destroyed. 

To be more precise, suppo se we have a dynamical sce- 
nario as described in Section II A and the 'out' dynamics 
A° has a basin such that the largest neighbourhood of 
A" contained in the basin has radius p > 0, and the dy- 
namics is uniformly contracting onto A° in the tangential 
direction at a rate r\ < 0. We can model the approach to 
A° along its weak stable manifold by a map x n +i — e v x n 
where x corresponds to the distance from A°. Perturbing 
this map by i.i.d. noise £ n that is uniformly distributed 
in [— er\/3, cr\/3] (such that the variance is a) we obtain 
an iterated function system of the form x n +\ — e v x n +£ n . 
We can see that fluctuations will drive x n to exceed p > 
if 



a > p(l-eV)/V3. 



(53) 



Consequently we expect that the 'out' phase and 'in' 
phase can no longer be distinguished once the noise has 
reached the order of this threshold. 



IV. NUMERICAL RESULTS AND SCALINGS 

In order to test our model of in-out intermittency 
(with and without noise) we consider a simple model 
mapping of the plane introduced in |10|j 

/Or, y) = Ml -x) + sxy 2 , 1.82e~ x y - y 3 ), (54) 

which has two parameters r € [0,4] and s S R. We 
can view this as a map of R 2 to itself that leaves Mi = 
R x {0} invariant. If s = 0, the map has the form of 
a skew product over the dynamics in x, i.e. it can be 
written as 



f(x,y) = (h(x),g(x,y)) 



(55) 



with 



h(x) = rx(l ~ x) and g(x, y) = 1.82e~ x y - y 3 , (56) 

where x € Mi . If we fix r and vary s we see that the latter 
two parameters do not affect the map restricted to Mi 
and so are normal parameters for the system restricted 
to Mj. 

An example of the behaviour of the transverse and tan- 
gential L.E.s around a window of periodicity for which 
the map (54) shows in-out intermittency is depicted in 
Figure |. We have also shown in Figure || a typical 
time series corresponding to the in-out intermittent be- 
haviour produced by this map. The top panel clearly 
demonstrates windows of periodic locking (correspond- 
ing to the 'out' phases), interspersed by chaotic windows 
(corresponding to the 'in' and 'reinjection' phases). One 
can also clearly see from the bottom panel the exponen- 
tial growth in the amplitude of the transverse variable y 
during the 'out' phases. 
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To study the effects of noise on in-out intermittency, 
it is informative to compare it with the analogous stud- 
ies of on-off. To do this we chose two sets of values for 
the control parameters r and s (namely r = 3.8800045, 
s = -0.3 and r = 3.82786, s = 0) in the map (||), 
corresponding to in-out and on-off intermittencies re- 
spectively. We perturbed the map with uniform noise on 
[0, Cl for the x dynamics and [— £, £] for the y dynamics. 
We have considered the two cases above where the noise 
is imposed (I) on the tangential variable x and (II) on 
the transverse variable y. 



Estimating the parameters for the noise-free 
model 



Non-normal parameter r 

FIG. 2: L.E. Ai in the tangential direction and At in the 
transverse direction for the map ( |54[ ) with s = —0.3 and vary- 
ing r. Initial conditions chosen to evolve to the attractor in 
Mi, i.e. with y — 0. Observe the 'periodic window' coin- 
cides with transverse instability of the attractor in Mi. In 
this region one can find in-out intermittent attractors for the 
full map that are not contained in Mi. Note that the at- 
tractor with transverse L.E. approximately —0.005 remains 
as a chaotic saddle during the periodic window; this saddle 
controls the in-phase. 
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FIG. 3: Example of time series generated by iterating (m) 
from a randomly chosen initial condition, after transient ef- 
fects have died out. The bottom plot shows y in a loga- 
rithmic scale, and the top second and the middle plots show 
x and y on linear scales. Parameter values are fixed at 
(r,s) = (3.88615,-0.3). 



Observe that the noise-free model for in-out inter- 
mittency, after a suitable non-dimensionalisation of the 
transverse variable has four parameters; A, 0, e, c with 
an arbitrary choice for application of the boundary con- 
dition at y = 1 after a suitable rescaling of y. We esti- 
mate these parameters as follows. We take a trajectory 
such that the transient has decayed and we can identify 
parts of the trajectory as 'in' or 'out' phase by choice of 
a suitable U°. For the map Q with r = 3.8800045 and 
s = —0.3, we can dchnc the trajectory (x n ,y n ) as being 
in the 'out' phase if 



Pi 



i = 1 



12} < 10" 



i.e. if it approaches the period 12 attractor {pi : i — 
1 . . . 12} for y = to within 10~ 4 . Using this criterion, 
we have depicted in Figure || the values of the transverse 
variable y at the entrance and exit of the 'out' phases 
identified by this procedure. Note that the exit point 
is more or less constant while the entrance point is dis- 
tributed in a manner consistent with exponential distri- 
bution in. z = — In. \y\ (see Figure |^). 



1. Estimating e. 

Since our model implies that probability decays 
from the 'in' phase at a rate e per unit time, this cor- 
responds to an exponential distribution of lengths of the 
'in' phases with average length of 'in' phase being 1/e. 
Hence 



1 

e = z? 

where one can easily approximate the quantity 
Li = {Average length of 'in' phase}. 



(57) 



(58) 



For the map @ with r = 3.8800045 and s = -0.3, we 
estimate L t ~ 5600 ± 200 and so e - 0.00018 ± 5 x 10~ 6 ; 
see Figure [| 
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FIG. 4: A plot of ln(y) at the entrance at the nth period, 
showing the value at the start of the 'out' phase (lower point) 
and that at the exit from the 'out' phase (upper point) for 
in-out intermittency, plotted against event number. Observe 
that exit never occurs closer than a certain distance from the 
invariant manifold at In \y\ — — oo. The events appear to be 
independent and uniformly distributed in time. 
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FIG. 5: A histogram of the entrance to the 'out' phases and 
the exit from the 'out' phases for in-out intermittency for the 
trajectory considered in Figure ^[ 



2. Estimating c. 

This is simply the largest transverse L.E. for A°, and 
can in some cases be obtained analytically. It can also 
be estimated numerically as the average growth rate of 
\y\ during the 'out' phases. For the example of in-out 



FIG. 6: The histogram of the lengths of the 'in' phases for the 
trajectory considered in Figure ||. This was used to calculate 
e from the inverse of the average length of the 'in' phases r. 



intermittency discussed above we can compute c to be 

0T434 =0Q1195±10 -b (59) 
12 v ' 

which is the transverse L.E. of the attracting period 12 
orbit in Mj. 



3. Estimating A and (3. 

These parameters can be obtained by examining the 
average rate of growth during 'in' phases. More precisely, 
we pick a threshold z t h = — hi(j/ih) which is large and 
time T > 0, and then examine all instances where the 
trajectory starts in the 'in' phase at z(in) = z t h and 
remains in the 'in' phase for at least a time T. 

Note that T needs to be chosen so that z(t) does not 
get too small (i.e. y{t) does not get too large) for to < t < 
to + T and one needs to be careful to avoid limiting the 
trajectory in such a way that may condition the mean 
or variance we are trying to measure, for example by 
choosing the threshold in y to be too large, or by choosing 
T to be so long that one will enter the nonlinear range. 

Subject to this, we can approximate A as the average 
value of 



(Hy(to + T))-Hy(t )))/T 



(60) 



over this ensemble of 'in' phases. Similarly (3 can be 
found as the standard deviation of this quantity from its 
mean value, per unit time. For the example of in-out 
intermittency in map ( |54j ) discussed above, we used up 
to 20,000,000 points of the trajectory, with T = 100 and 
an ensemble of in-phase segments of the same trajectory 
with z t h ~ 30 to find that 



A 



-0.0042, /3 ~ 0.0135. 



(61) 
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where there is an expected maximum error of approxi- 
mately 5% (see Figure 0) . 



B. Lack of fit to a Fokker-Planck model of on— off 
intermittency 
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FIG. 7: Histogram of the distribution of calculated values 
of A and estimation of A and (3 from the in phases for the 
trajectory considered in Figure ^| (see text for details). 



4- Check: an independent estimate of y,. 

Recall that the ratio of the times spent in the 'in' to 
the 'out' phases can be obtained from the stationary dis- 
tributions in the form 



Rio = 



<j>o 



.4 
B 



Cft 

e 



Using our knowledge of e and c we can easily obtain fi 
and check this against the theoretical prediction (20), 
given approximations of the quantities $' = { Asymp- 
totic proportion of time spent on the 'in' phase } and 
$° = {Asymptotic proportion of time spent on the 'out' 
phase} = 1 — P l . For the case of in-out intermittency 
considered above, the above estimates of the parameters 
imply that 

H 0.0451, 

which allow us to obtain a numerical estimate for Ri Q ~ 
3.02 with an estimated maximum error of 15%. 

Using instead the measured ratio of average length of 
'in' to 'out' phases we obtain an estimate for Ri D — 2.56 
with an estimated maximum error of approximately 10%. 
These two estimates of Ri a clearly agree to within esti- 
mated maximum error. 

Note that the neighbourhood U° of the 'out' dynamics 
A must be chosen such that the U° D Mj is forward 
invariant. It can be chosen as small as desired, though 
if it is very small then we will not recognize 'out' phases 
unless they come very close to Mj. 



It is interesting to note that it is not possible to fit 
in-out intermittent dynamical data by an on-off model. 
This is because the on-off model requires only two pa- 
rameters, the transverse L.E. A > and its variance per 
unit time /3 2 . If we examine the attractor in Mj, we can 
compute a positive L.E. (c above), but the variance would 
be zero. Alternatively we can compute the 'in' phases as 
discussed above and obtain both a A and a /3 2 > 0, but 
in that case A < 0. Thus either choice will be invalid. 

Alternatively one could compute A//3 2 from the scal- 
ing of the probability density near y = 0, but then it is 
not clear how to make a sensible choice for A or /3 and 
therefore we can determine only one of the parameters in 
the model. 



C. Probability distribution for the case with noise 

Figure || shows the results of three distinct ways of cal- 
culating the asymptotic probability distribution function 
(pdf) of q(y); resp. s(y), for a fixed noise level (in this 



case (^y 



io- 



from the direct integration of (plf) us- 



ing the estimate values of the parameters e, c, A and /3; 
from the full analytical solution ( |3l| ) and finally from the 
direct numerical measure of the pdf of Q(y). 
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FIG. 8: Histogram of the pdfs q(y) on the in-phase and 
s(y) on the out-phase as a function of the transverse variable 
y for transverse added noise o~ v — 10 -8 for the examples of 
in-out intermittency in ( |54| ) discussed in the text. Note that 
the analytical solution joins the two regions, \y\ < a and 
\y\ > a I 'f3 by means of an internal boundary layer. Note that 
any algebraic scaling of q(y) m y a implies the same scaling 
for s(y). 

We have also plotted in Figure || the influence of differ- 
ent values of the noise level £ on the transverse variable 
y in the pdf of Q(y), for both in-out and on-off cases, 
and discuss the behaviour in the figure caption. 
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FIG. 9: Histogram of the pdf Q(y) as a function of the trans- 
verse variable y for transverse added noise for the examples of 
on-off and in-out intermittency in ([si]) discussed in the text. 
Observe that there is a power law scaling across a wide range 
of y that changes to a constant smooth density (plateau) for 
very small y. The turning point corresponds, as expected, 
to the noise level imposed; the response of on-off and in-out 
intermittency to the addition of transverse noise can be seen 
to be similar. 



D. Length of the average laminar phase as a 
function of noise 



FIG. 10: The difference between on-off and in-out intermit- 
tency shows itself qualitatively in the behaviour of the average 
laminar phase size as a function of the noise level £ = \J~Zo for 
noise addition to the tangential variable x. This illustrates the 
case for tangential and transverse added noise to map (^ij). 
The noise on x is uniformly distributed in the range [0,^], 
while that on y is in the range [— f y , +Cy]- The parameter val- 
ues are r = 3.8800045, and s = —0.3. The decay for noise on 
x in the in-out case corresponds to the ^-perturbation getting 
large enough to destroy the basin of attraction of the 'out' dy- 
namics; by contrast the addition of tangential noise to on-off 
intermittency only makes a difference at much larger levels. 



Noise on tangential variable 



We calculated for the map (54) the scaling of the length 



of the average laminar phases, as a function of the noise 
level (,. Our results are depicted in Figure [l(]. As can be 
seen, in this case there is a significant difference between 
the on-off and in-out cases. While both show very simi- 
lar behaviours for low noise levels, at higher noise levels, 
the behaviours show some distinct differences. In partic- 
ular, average laminar sizes corresponding to the in-out 
case drop off rather suddenly, whereas for the on-off case 
there is an increase in the size of the average laminar 
phases before it too drops off sudd enly. This corresponds 
well with the discussion in Section III B , where we argued 



that in the case of in-out, a sudden change in dynam- 
ics would be expected at a certain noise level. For the 
map considered here, we observe that the local basin of 
the 'out' dynamics is of the order 10 -6 which from ( |5^ ) 
would suggest a noise threshold that is at most 10" 6 . 

Also shown on this Figure are the results for the case of 
on-off intermittency for the same map at a different pa- 
rameter value (see caption) . One would expect a decrease 
in the average length of laminar phases until a is of the 
same order of the threshold defining the laminar phases, 
which in this case is ~ 10~ 3 . This is to be contrasted 
with the in-out case, where the drop off is more sudden 



and occurs at much lower noise levels (i.e ~ 10 -6 ). This 
level of noise seems to be enough to disrupt the periodic 
attractor in the invariant submanifold and it is interest- 
ing to note that the noise level at which this occurs is of 
the same order of magnitude as the size of the parameter 
window of periodicity. 



2. Noise on transverse variable 



For the noise on the transverse variable, on-off and in- 
out behave quite similarly in terms of the average laminar 
phases, showing a smoother decay than the case (I), with 
the dramatic drop occurring around 10~ 4 , which is closer 
to the threshold level. As can be seen from Figure O, the 
'out' chains, initially dominant, decay rapidly, while the 
'in' chains are on average of the same size for a wide band 
of noise levels a, up to values of about 10~ 6 . Note also 
that the actual percentage of time spent in the 'in' chains 
actually increases at high noise levels before decaying to 
zero, while the time spent in the 'out' chains decreases 
monotonically. 
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FIG. 11: Comparison of the average length of laminar phase, 
'in' phase and 'out' phase as a function of the noise level a 
for addition of transverse noise in the variable y; parameter 
values are r = 3.8800045 and s = -0.3. 



3. Average mean crossing time through y = 

One can similarly analyze for the model the average 
mean crossing time M through y — in presence of noise 
on y (case (II)). For positive transverse L.E. (At > 0), 
that is for parameter values away from the blowout point, 
one expects for the case of on-off a typical growth given 
by f§ 

M » £— a' 6 + O(l). (63) 
A 

For the on-off case, |l8) predicts that at blowout point 
the scaling has the asymptotic form (A — > 0) 

M«-^+0(l). (64) 

In the in-out case we find the scaling for small \a\, can 
be approximated from the expansion of j5C| ) in a, con- 
sidering At ^ 0, in the form 

M fa Ci + C 2 In a + C 3 In 2 a + C>(ln 3 a) . (65) 

We verified this for the case of in-out in the map above, 
by using the same parameter values, except for the nor- 
mal parameter b which was chosen such that At = 
ln(1.82) + b(x) r = 0, in order to enable us to calculate 
numerically the dependence of the average mean crossing 
time M through y = for the blowout point (At = 0) . In 
Figure [jj] we verify the above predictions of the scaling 
for the mean crossing time M for the three cases At < 0, 
At = and At > 0. Note in particular the case At > 0, 
which in contrast with the on-off case needs the term in 
In 2 a to be included (see Ashwin and Stone |Q for the 
on-off version of Figure fi^) . 



FIG. 12: Mean crossing times of variable y through zero, as 
a function of noise level a, with noise added on the transverse 
variable y, for three cases: At < 0, At = 0, At > 0. Parame- 
ter values are r = 3.8800045, v = 1.82, a = -1 and s = -0.3. 
The parameter b is varied to obtain normal variations on At- 

4- Probability distributions of laminar phases 

We have as yet been unable to compute a closed form 
approximation for the probability distribution of the 
laminar phases for in-out intermittency for the Fokker- 
Planck model, but Figure |l3| suggests that the scalings 
are analogous to those obtained in |l(| for the discrete 
Markov model. Note the presence of an inflection point 
and 'shoulder' in the in-out distribution corresponding 
to a relatively high number of long laminar phases. This 
shoulder appears to persist on addition of noise. By con- 
trast, the distribution for on-off laminar phases does not 
show such a shoulder. 



V. DISCUSSION 

We have proposed a continuum model of the statis- 
tics of the transverse variable for in-out intermittency 
in the form of a Fokker-Planck model with delay integral 
boundary conditions to model the deterministic propaga- 
tion of probability density near the unstable manifold of 
the 'out' phase. This presupposes that the 'out' dynam- 
ics are periodic in the invariant manifold, but if they are 
not then similar models could be derived in the form of 
coupled Fokker-Planck equations. Such models are then 
well-adapted to model the addition of further noise. 

Although in-out intermittency has a number of simi- 
larities to on-off, we see that there are differences in their 
statistical properties. In particular, the addition of noise 
'tangential' to the dynamics can lead, at least in our ex- 
amples, to significant changes in behaviour at much lower 
noise levels than for on-off. 
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FIG. 13: Scalings of the probability distribution of laminar 
phases for (a) on-off and (b) in-out intermittency with noise 
added in the transverse direction. The noise perturbations on 
y are uniform in [— a y , +a y ]. Parameter values are (a) (r, s) = 
(3.82786,0) and (b) (r,s) = (3.8800045, -0.3). Observe that 
the on-off statistics limit do not show the presence of the 
inflexion point clearly visible in the in-out statistics. 



We have demonstrated how, given an in-out intermit- 
tent signal, it is not possible to sensibly fit the parameters 
to on-off intermittency from the dynamical data avail- 
able. There is clearly a lot more one could examine in 
such models, for example, the scalings of the variance and 
mean first crossing times with the various model param- 
eters and noise level; there is work presently in progress 
that aims to understand the variation of such scalings on 
change of system parameters. 



Implicit in our work here is the assumption that the 
in-out intermittent attractor supports a natural ergodic 
invariant measure, such that almost all points attracted 
to the attractor will display the same stationary statis- 
tical behaviour. Although we do not doubt this for the 
models considered so far, it does seem possible that in- 
out intermittency may give rise in certain circumstances 
to behaviour that is not ergodic, and one needs to bear 
in mind the possible existence of such behaviour. 
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